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We report experiments on the deformation and transport of an elastic fiber in a viscous cellular 
flow, namely a lattice of counter-rotative vortices. We show that the fiber can buckle when ap- 
proaching a stagnation point. By tuning either the flow or fiber properties, we measure the onset 
of this buckling instability. The buckling threshold is determined by the relative intensity of vis- 
cous and elastic forces, the elasto- viscous number S p . Moreover we show that flexible fibers escape 
faster from a vortex (formed by closed streamlines) compared to rigid fibers. As a consequence, the 
deformation of the fiber changes its transport properties in the cellular flow. 

PACS numbers: 



The interaction of a deformable body with a vis- 
cous flow is found in a wide range of situations, rank- 
ing from biology to polymer science. Cells or micro- 
organisms in water use deformation of high-aspect ratio 
flexible organelles called cilia or flagella to move at low 
Reynolds number Numerous studies have been ded- 
icated to locomotion and hydrodynamic interactions at 
low Reynolds number, including the study of biological 
0, Q or model systems , but a large number of ques- 
tions remains still to be answered. The stretching of a 
flexible polymer in a viscous flow is well studied @ . The 
so called coil-stretch transition is responsible for the spec- 
tacular macroscopic properties of polymer solutions as for 
example normal stresses or high elongational viscosities. 
This transition might have an analog for elastic fibers 
that buckle by viscous forces. Microscopic buckling of 
fibers has indeed been shown theoretically to lead to orig- 
inal rheological properties of their suspensions, such as 
the appearance of normal stress differences P3]. Another 
fundamental question is the modification of the transla- 
tional dynamics of an object induced by its deformation 
in a viscous flow as studied by Young and Shelley [TTj . 

The complex dynamics of such systems originates from 
different sources: first, the dynamics is dependent on the 
size of the deformable body with respect to the charac- 
teristic lengthscalcs of the flow. Second, when the struc- 
ture is deformed by viscous forces, the local surrounding 
flow will be modified leading to non-linear effects. Hy- 
drodynamic interactions between flexible objects further 
increase complexity [l2j . The experimental study of sim- 
ple model systems, such as homogeneous clastic fibers, 
evolving in a controlled flow geometry, could help under- 
standing these complex interactions. 

In this communication we experimentally study the de- 
formation of an isolated macroscopic filament approach- 
ing a stagnation point and the subsequent modification of 
its transport through a cellular flow of counter-rotating 
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FIG. 1: Successive shapes of two fibers experiencing the same 
flow field but with different Young moduli (At = OAs, from 
red to purple). The more rigid fiber (bottom left) does not 
buckle whereas buckling is observed for the more flexible fiber 
(up right). The background flow is visualized independently. 



vortices (see figure [T]). First we characterize the buckling 
transition of the fiber varying the properties of the fiber 
and the viscous flow independently and identify the buck- 
ling threshold. The threshold is shown to be function of 
the relative intensity of the viscous and elastic forces, the 
so called S p number We then discuss differences in 
the transport of rigid and flexible fibers in the cellular 
flow and show that deformation significantly changes the 
transport properties. 

The lattice of counter-rotating vortices is obtained us- 
ing electromagnetic forcing [f4|, Ell ■ An electric current is 
applied to an electrolyte placed in a rectangular cell (Fig. 
(2^). A lattice of 4x5 magnets (NdFeB, Supermagnets) 
with alternating orientation is placed underneath the cell, 
creating a spatially periodic magnetic field. Locally vary- 
ing Lorentz forces act on the flowing electrolyte leading to 
an array of counter-rotating vortices of size W = 3 cm. 
The electrolyte is a 50:50 mixture of Polyethylenegly- 
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FIG. 2: a) Experimental setup: a lattice of magnet with al- 
ternated orientation is placed below an electrolyte, b) Veloc- 
ity field. Closed circles represent stagnation points; crosses, 
centers of the vortices, c) Maximum velocity Uo vs electric 
current /. The line is the best linear fit. d) Color mapping 
of the local compression rate e in a cell. The dashed line is 
one experimental trajectory of the fiber. The dotted circle 
represents 0.8 e max . 



col 1000 (Fluka) and purified water, into which NaCl is 
added up to saturation (viscosity, 77 = 35 ± 5 mPa.s, sur- 
face tension 7 = 43 ± 1 mN/m. The electric current 
/ = —400/ + 400 mA is applied via a DC power supply 
using carbon fiber electrodes (Toyobo). 

The array of vortices generates a cellular flow and an 
array of stagnation points. Close to the stagnation points 
the flow can be described by hyperbolic streamlines, cor- 
responding to a purely elongational flow with constant 
elongation/compression rate e. Close to the center of 
the vortices the flow can be described by solid rotation 
without any deformation of the fluid. Note that the flow 
consists of closed streamlines and small particles can not 
be transported across the flow field when diffusion and 
particle inertia are neglected. 

The velocity field is determined by Particle Imaging 
Velocimetry (using DAVIS software) and presented on 
figure [2j All our experiments are performed in a selected 
region, a 3x3 lattice of stagnation points, where boundary 
effects are negligible. The velocity field is well described 
by the following equation: 

/ sin(x)cos(y) \ , 1N 
U(x,y) = U V\ • r~\ (1) 

v ,yj \ -cos(x)svn(y) J v ' 

where i (resp. y) is the dimensionless coordinate irx/W 
(resp. ny/W). x and y are set to zero at a stagnation 
point. The parameter Uq is found to be proportional to 
the electric current / (see fig. &p). Note that in the range 
of velocities used the Reynolds number is always below 



5. The local compression rate of the flow reads: 

V, I) = ± nU ^ cos(x)cos(y) (2) 

with opposite signs considering either compressive or ex- 
tensive quadrants. The direction of compression remains 
within one quadrant always the same. The compression 
rate is locally variable, being maximum at the stagna- 
tion point (e max = ttUq/W), and zero at the center of 
the vortices. 

During its motion in the lattice, a fiber thus ex- 
periences a variable compression rate e(x,y). This 
is illustrated in figure [2}i, where the local compres- 
sion rate is color mapped in a cell centered on a 
stagnation point. All our experimental trajectories 
reach 80% of e max . Thus, for sake of simplicity, we will 
hereafter use this value to characterize the flow influence. 

Experiments are performed with macroscopic fibers 
(length L=12mm ± 0.5 mm, radius r=85 ± 5 /im) made 
by crosslinking vinylpolysiloxane with a curing agent 
(Zhermack, 8 Shore A) in a glass capillary. The filament 
is then carefully extruded from the capillary. Optical 
measurements on a microscope have been used to mea- 
sure the radius of the fiber as well as to verify that it is ho- 
mogeneous and not damaged over the whole length. The 
ratio between its length and the cell size a = L/W = 0.4 
is fixed for all experiments. The mass ratio between the 
curing agent and the polymer can be varied from 1:2 to 
3:1. The Young's modulus is obtained via rhcological 
measurements using Y — 3G' and can be tuned from 
y=75 kPa to Y = 180 kPa. 

The density of the fiber is lower than the density of the 
fluid (Ap = 0.15.10 3 kg/m" 3 ) and the fiber thus floats at 
the surface. More precisely, the equilibrium position of 
the fiber at the interface is determined by the interplay 
between gravity and surface tension [l6[ . For our system, 
surface tension dominates, leading to the absence of a 
meniscus. The degree of immersion of the fiber is than 
solely controlled by the wetting angle between the liquid 
and the fiber. We find this angle to be 90 ± 5 deg and 
the fiber is approximativcly half-immersed. As a conse- 
quence, the motion as well as the deformation of the fiber 
are bidimensional, constrained at the surface of the flow. 
As the interface is not deformed, the influence of surface 
tension on the dynamics of the filament is negligible. 

Pictures of the fiber are taken with a digital camera 
(PixelLink, 1024x768, 10 fps). Using standard detection 
procedures (Image J software), we measure the position 
of the center of mass of the fiber. A fit of the shape is 
performed with spline functions (Numerical Recipes 3.3), 
and used to extract the mean curvature of the fiber as a 
function of time t, C(t). 

We will now use the experimental set-up to study the 
buckling instability of the flexible fibers in the vicinity of 
a stagnation point. Figure Q] shows two typical trajecto- 
ries of a fiber with different elastic moduli, but evolving 
in the same flow. One observes that the rigid fiber does 



3 



not buckle when approaching the stagnation point, the 
flexible fiber however is deformed by the flow and shows 
a buckling instability. 

To characterize this buckling instability in detail we 
construct the probability of a fiber to buckle in a given 
flow. We consider trajectories of the fiber cell by cell, 
where a cell is a square of width W, centered on a stag- 
nation point. One trajectory is defined as the path in 
the vicinity of one stagnation point as illustrated by the 
dashed line on figure [5Ji. The fiber is said in "coiled 
state" if the mean curvature is larger than a given thresh- 
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0.125 mm , i.e. to a radius of curvature 



of approximatively 2L/3. This corresponds to almost a 
quarter of a circle. No quantitative modifications of the 
results are observed if the threshold is slightly tuned. We 
do not consider trajectories where the filament enters the 
cell in a coiled state. We determined, for each experiment 
the ratio between the number of trajectories displaying a 
coiled state N C oiied and the total number of trajectories 
Ntot- Assembling over 20 experiments, around 7000 tra- 
jectories of the fiber are considered, with a total number 
of buckling events equal to «500. 

The fiber is predicted to buckle when viscous forces 
overcome elastic forces. Their relative intensity is given 
by the Sperm number S p [HI]: 



(3) 



Our set-up allows to vary the properties of the fiber 
and of the viscous flow separately tuning the two forces 
independently. The inset of figure[3]shows the probability 
to buckle as a function of the maximum velocity of the 
flow U for different Young's moduli (Y=75, 120 and 180 
kPa). One observes for all cases an increase with Uq. 
One observes furthermore that the probability to buckle 
is higher for the flexible fibers. 

Figure [3] now shows the probability to buckle for all ex- 
periments and its average as a function of S p . Within the 
experimental resolution the data points fall on a single 
curve, showing that the S p number is indeed the control 
par ameter of the system, as predicted by Young et al. 

The probability to buckle is zero below S p = 120, it 
increases at S p ~ 120 — 150 and reaches at our highest S p 
about 20%. Theoretical predictions from linear stability 
analysis in a hyperbolic flow show a threshold of S* =153 
for the buckling instability 17j [ll| . Note that their defi- 
nition of S* includes the logarithmic correction from slen- 
der body theory and reads S* = Y (in(L < /R)-i/2) 7*- We 
correct the experimentally obtained value for the buck- 
ling threshold in the same way and find S* ~400. In 
the experiments the fibers are only half-immersed and 
viscous forces are certainly smaller then those used in 
the theoretical analysis. This is likely to be the main 
reason for the higher threshold observed in our experi- 
ments. Furthermore, even if we only consider filaments 
passing close to the stagnation point, the filament is of 



significant length compared to the typical lengthscale W 
of the cellular flow and the approximation of hyperbolic 
streamlines docs not hold in our experiments. 

We do not observe an abrupt increase of the probability 
towards 100% at the threshold but a continuous increase 
with S p . At the highest value of S p accessible in our 
experiments the filament still only buckles one out of four 
times it passes a stagnation point. This indicates that 
even if the threshold of buckling is controlled by S p the 
dynamics in the cellular flow are more complicated. It 
is likely that the probability might increase further when 
going to higher values of S p . 
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FIG. 3: Probability of buckling vs S p for different Young 
moduli (red cross Y=75kPa, blue triangle Y=120 kPa, green 
open circle Y=180 kPa). Black closed circles represent the 
mean value. Inset: Probability of buckling as a function of 
U . 



With our experiment we can also address an other fun- 
damental question which stands in the modification of 
the translational dynamics of the filament induced by its 
deformation. To do so we study the escape of a filament 
from a vortex. The finite size of the fiber already changes 
its dynamics in the cellular flow compared to a pointlikc 
particle (see figured}. We will thus compare deformable 
fibers (S p > 120) to rigid fibers (S p < 120). The fibers 
are placed in a vortex parallel to the axis of compression 
and at a given distance b (see figure [4jj). We then record 
the number of trajectories the filament needs to escape 
from the vortex and the number of buckling events during 
this escape for over 250 experiments. 

Figure [4] shows the number of trajectories as well as 
their average above and below the buckling threshold 
S p = 120 as a function of b normalized by the size of 
the vortex W/2. Note that the range of b studied cor- 
responds to trajectories of the filament reaching at least 
80 % of the maximum compression rate, identical to the 
experiments used to study the buckling probability. 

Trivially, one observes that the filament needs longer 
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FIG. 4: Number of trajectories before escape as a function 
of b/(W/2). Open symbols, individual experiments; closed 
symbols, the mean value. Red circles S p > 120, blue squares 
S p < 120. Inset, definition of b. 



to escape from a vortex if it is further away from the axis 
of compression. One observes also that even the rigid 
filaments escape from the vortices after a finite number 
of trajectories, showing the importance of the effect of 
the finite size of the filament. The number of trajectories 
needed to escape from a vortex fluctuates strongly in the 
case of the rigid fibres, indicating that small changes in 
the initial position (not measured in this analysis) are im- 
portant for the dynamics. Less fluctuations are observed 
for the flexible filaments. When comparing the average 
number of trajectories of the flexible fibers and the rigid 
fibers one observes that the flexible filaments escape sig- 
nificantly faster. This is a clear signature of the effect of 
the filament's deformation on its transport properties. 
Young and Shelley show that the onset of the buck- 



ling instability is coupled to the onset of a translational 
" Brownian like" motion of the fiber across the lattice of 
stagnation points [ll|. Above the stretch coil transition, 
fibers present a diffusive behavior, with however, a non 
trivial dependence of the diffusion coefficient with the rel- 
ative intensity of viscous and elastic forces. The motion 
of a fiber in one cell can be seen as a 'collision' with the 
stagnation point. Sequences of these collisions can then 
lead to diffusion of a fiber. Our results are the first indi- 
cation that flexible fibers indeed diffuse faster than rigid 
ones through the cellular flow. 

In conclusion, we have built an experimental setup al- 
lowing to study the deformation and the transport of a 
ccntimctric fiber in a cellular flow, formed by counter- 
rotating vortices. We have shown that the elastic fiber 
can buckle due to the compressive viscous forces. We 
have identified the threshold for the buckling instability, 
which is function of the so called S p number, being a 
measure of the relative intensity of viscous and clastic 
forces. These experimental observations are in reason- 
able agreement with linear stability analysis in a purely 
hyperbolic flow. We have also studied the translational 
dynamics of the filament in the cellular flow field. We 
have shown experimentally that fibers can escape from a 
given vortex due to their finite size compared to the size 
of the vortex. Flexible fibers are shown to escape after 
significantly less tours compared to rigid fibers. This is a 
first evidence that the deformation of a filament changes 
its transport properties in the cellular flow. 
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